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Ionospheric  Modeling 


I.  INTRODUCTION 

The  reflection,  transmission,  refraction,  and  scattering  characteristics  of  the 
ionosphere  have  been  used  for  radio  communication  since  the  discovery  of  the  ionos¬ 
phere  in  the  earth's  atmosphere.  For  successful  radio  communication,  it  is  essen¬ 
tial  to  predict  the  behavior  of  the  ionospheric  region  that  will  affect  a  given  radio 
communication  circuit.  Such  a  prediction  will  identify  the  time  periods,  the  path 
regions  and  the  sections  of  high  frequency  bands  which  will  allow  or  disrupt  the  use 
of  the  selected  high  frequency  communication  circuit.  This  need  for  prediction 
leads  to  modeling  of  the  ionosphere. 

A  model  is  a  numerical  statistical  description  of  the  ionosphere  in  terms  of 
location  (geographic  or  geomagnetic  latitudes  anc  longitudes,  time  (solar  zenith 
angle),  seasons,  and  other  factors  such  as  the  solar  activity  (10.  7  cm  flux,  sunspot 
numbers).  The  empirical  equations  are  derived  from  the  dependence  of  the  observed 
phenomena  on  various  parameters  mentioned  above.  These  observed  phenomena 
include:  the  behavior  of  critical  frequencies  foE,  foF^,  foFg*  and  foEs  for  the  E, 
F.,  F„,  and  sporadic  E  layers;  the  altitudes  for  peak  (maximum)  electron  densities 

(Received  for  publication  26  January  1982) 

The  critical  frequency  is  the  limiting  radio  frequency  below  which  a  radio  wave 
is  reflected  by,  and  above  which  it  penetrates  and  passes  through,  the  ionized 
medium  (an  ionospheric  layer)  at  vertical  incidence. 
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for  the  layers  h  E,  h  i-'  ,  and  h  F,  and  the  half-thickness  widths  y  E,  v  1',; 
and  y  These  models  are  called  phenomenological  models. 

The  IXS-78  model,  based  on  the  analysis  by  Jones  et  al6  of  world- wide, 
ground-based  ionosonde  data,  predicts  only  the  bottomside  of  the  ionosphere.  The 
Bent  model  predicts  the  total  electron  content  of  the  ionosphere  in  the  altitude 

range  from  150  to  2000  km,  without  a  direct  consideration  of  the  bottomside  E  and 

3 

Fj  layers.  The  Ching  and  Chiu  model  covers  the  altitude  range  from  110  to 
1000  km.  Instead  of  parabolic  layers1  they  assume  Chapman  functions  for  the  elec¬ 
tron  density  distributions  in  the  E,  F . ,  and  F0  layers.  Later  Chiu4  modified  the 
3  l  £ 

Ching-Chiu  model  to  incorporate  the  polar  ionosphere.  Their  models  are  useful 
only  for  studying  the  large  scale  phenomena  such  as  global  thermospheric  and 
ionospheric  calculations. 

Using  ionospheric  data  from  ESRO  satellites,  Kohnlein5  extended  the  altitude 
range  up  to  3500  km.  He  suggested  a  "differential  approach"  for  the  ionos¬ 
pheric  modeling.  He  separated  small  scale  special  structures  such  as  the 
equatorial  trough,  the  mid-latitude  trough  and  the  polar  ionosphere,  from  the  large 
scale  global  structure.  He  modeled  these  individual  structures  and  added  them  into 
the  global  structure.  His  method  reduces  the  number  of  coefficients  otherwise 
needed  to  model  the  complicated  ionospheric  behavior. 

The  other  approach  for  ionospheric  predictions  is  to  use  theoretical  models.  8-12 
These  are  based  on  the  physical  processes  responsible  for  the  production,  main¬ 
tenance  and  decay  of  the  ionosphere.  A  theoretical  model  would  thus  rely  on  the 
process  of  ionization  of  neutral  atmospheric  constituents  by  the  incident  solar  ex¬ 
treme  ultraviolet  radiation,  the  transport  processes  such  as  diffusion  and  neutral 
winds,  and  also  on  the  effect  of  electric  and  magnetic  fields  on  the  transport  pro¬ 
cesses.  Essentially  the  theoretical  model  tries  to  explain  the  experimental  observa¬ 
tions  in  terms  of  known  physical  processes.  In  addition,  this  approach  seeks  new 
physical  processes  to  explain  the  differences  between  the  observational  results  and 
the  predictions  based  on  the  theoretical  considerations. 

Radio  communication  can  be  divided  in  two  main  categories.  Ground-to-ground 
radio  communication  is  based  on  the  reflection  and  scattering  characteristics  of  the 
ionospheric  layers.  On  the  other  hand,  ground  to  satellite,  satellite  to  ground,  or 
satellite  to  satellite  radio  communications  rely  on  the  transmission  and  refraction 
characteristics  of  the  ionosphere.  The  main  goal  of  any  modeling  effort  is  to  pre¬ 
dict  the  periods  of  good  or  poor  radio  communications  for  the  selected  paths  to  en¬ 
able  a  continuous  undisrupted  communication  through  the  ionosphere  or  by  some  other 
means. 

The  half-thickness  width  ym  of  the  ionospheric  layer  is  determined  with  the 
assumption  that  the  layer  has  a  parabolic  shape.  * 

(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 

See  References,  page  25. ) 
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In  the  following  we  will  consider  several  ionospheric  models  which  are  routine¬ 
ly  used  (or  are  available)  for  the  prediction  and  specification  of  the  ionosphere. 

The  emphasis  here  is  on  acquainting  the  user  with  the  modeling  programs  and  their 
limitations.  This  report  in  no  way  attempts  to  review  the  scientific  literature  for 
a  determination  of  the  state  of  the  art  of  modeling  efforts.  Therefore  only  the 
essential  references  will  be  cited  in  the  following  text. 

First,  we  will  consider  the  numerical -phenomenological  models.  Then  we  will 
consider  the  theoretical  models.  This  will  be  followed  by  the  modifications  to 
models  to  take  into  account  high  latitude  phenomena  such  as  the  auroral  E  layer  and 
the  mid-latitude  F  region  trough.  In  the  conclusion  section,  we  will  look  at  the 
limitations  of  these  models  and  a  possible  approach  to  overcome  these  limitations. 


2.  THE  NUMERICAUPHENOMENOLOGICAL  MODELS 

At  present  the  three  most  widely  used  numerical  models  for  ionospheric  pre¬ 
dictions  are:  (1)  The  ITS-7  8  model,  (2)  The  Bent  model,  and  (3)  The  Ionospheric 
Communications  Analysis  and  Prediction  Program  (IONCAP).  In  addition,  the 
4-D  model  of  the  Air  Force  Global  Weather  Central  and  the  Bradley  model  will  be 
considered. 

2.1  The  ITS-78  Model 

The  main  purpose  of  this  model  *  is  to  predict  long  term  performance  of  com¬ 
munication  systems  operating  in  the  2-30  MHz  frequency  range.  The  model  serves 
an  essential  function  of  assistance  in  the  planning,  design,  and  operation  of  high 
frequency  (HP")  communication  systems. 

The  model  provides  (output)  circuit  operational  parameters  such  as  the  maxi¬ 
mum  usable  frequency  (MUF),  optimum  traffic  frequency  (FOT),  and  the  lowest 
usable  frequency  (LUF).  In  addition  to  the  regular  E-layer  propagation  mode,  it 
takes  into  account  propagation  via  the  sporadic  E  layer.  The  program  computes 
all  the  probable  modes.  It  computes  the  system  performance.  For  that  purpose  it 
calculates  the  antenna  patterns  and  gains  for  10  most  commonly  used  antennas.  It 
also  has  a  program  to  determine  MUF  as  a  function  of  the  magnetic  activity  index  Kp. 

The  inputs  for  the  ITS-78  model  are  the  date.  Universal  Time,  geographical 
location  (latitude  and  longitude)  of  the  transmitter  and  receiver,  and  sunspot  number. 
To  compute  the  system  performance  the  model  needs  the  antenna  parameters,  the 
radiation  power  of  the  transmitter,  and  the  signal  to  noise  ratio  of  the  receiver. 

The  ITS-78  model  and  its  computer  program  was  developed  by  the  Institute  of 
Telecommunication  Sciences,  ESSA,  Boulder,  Colorado.  The  model  is  based  on 
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the  presentation  of  the  ionospheric  characteristics  in  a  form  of  synoptic  numerical 

13  6 

coefficients  developed  by  Jones  and  Gallet  and  improved  by  Jones  et  al.  The 
important  features  of  the  ITS-78  model  are  the  parameters  for  the  D,  E,  Eg 
(sporadic),  and  F  ^  layers  of  the  ionosphere. 

For  the  D  region  the  model  considers  only  the  absorption  losses.  The  non- 
deviative  absorption  is  in  the  form  of  a  semi -empirical  expression.  It  enables  the 
user  to  compute  the  losses  for  the  HF  frequencies  penetrating  the  D  layer.  The 
deviative  absorption  losses  are  included  in  the  loss  calculations  as  uncertainty 
factors. 

For  the  E  region  the  model  computes  the  parameter  foE.  It  assumes  a  con¬ 
stant  height  of  110  km  for  the  maximum  (peak)  electron  density  of  the  E  layer,  with 
a  constant  semi -thickness  of  20  km.  The  numerical  coefficients  for  foE  are  based 
on  the  experimental  ground  ionosonde  data  during  high  solar  activity  phase  in  1958, 
and  the  low  solar  activity  phase  in  1964. 

For  the  F^  region  the  model  computes  the  parameters  foFl?,  the  height  of  maxi¬ 
mum  electron  density  h  F  and  the  semi -thickness  y  F0  of  the  F  layer.  These  are 
in  the  form  of  numerical  coefficients  for  the  high  (1958)  and  low  (1964)  phases  of 
solar  activity.  Both  the  E  and  F^  layers  are  assumed  to  be  parabolic  in  shape. 

The  sporadic  E  (Eg)  layer  could  be  very  helpful  or  harmful  to  radio  communica¬ 
tions  depending  on  the  nature  of  the  E  layer.  A  blanketing,  totally  reflecting  E 

s  s 

laver  extends  the  frequency  range  of  the  E-mode  communications.  However  a 

semi-transparent  or  partly  reflecting  layer  would  cause  serious  multipath  and 
mode  interference  and  would  be  detrimental  to  communication  systems.  Using 
numerical  coefficients,  the  ITS -7  8  model  computes  foE,  only  for  the  ordinary  ray. 
(The  earth's  magnetic  field  splits  the  incident  ray  into  the  ordinary  and  the  extra¬ 
ordinary  rays.  )  The  numerical  coefficients  are  for  both  the  high  (1958)  and  low  (19  64) 
phases  of  solar  activity.  As  the  model  does  not  predict  the  occurrence  of  E  ,  the 
foEg  maps  are  used  only  when  propagation  via  regular  E  layer  is  not  possible.  To 
compute  the  system  performance  the  model  incorporates  three  kinds  of  noise: 
galactic,  atmospheric,  and  man-made. 

To  determine  the  operational  parameters  such  as  the  maximum  usable  frequency 

(MUF),  the  model  computes  the  path  geometry  (between  the  transmitter  and  the 

receiver).  The  parameters  computed  in  the  path  geometry  are  the  path  length,  path 

bearing  (azimuth),  and  the  solar  zenith  angle  \  of  the  sun.  The  model  computes 

paths  for  reflections  from  the  E,  E  ,  and  layers.  These  are  (  ailed  the  E.  E  . 

s  2  s 

and  F2  modes.  The  paths  could  involve  one  layer  (single  mode)  and  one  reflection 
(one  hop),  or  would  involve  more  layers  (multiple  modes)  and  more  reflections 
(multihop). 


13.  Jones,  W.  B.  ,  and  Gallet,  If.  M.  (1960)  Ionospheric  mapping  by  numerical 

methods,  TTU  Telecomm.  J.  ,  12:260. 

- -  ^ 
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For  considering  the  wave  propagation  the  electron  density  distribution  with 
altitude  is  needed.  Both  the  F  and  F  layers  are  assumed  to  have  parabolic  shapes. 
The  maximum  usable  frequency  (MUF)  is  obtained  by  multiplying  the  critical  fre¬ 
quency  of  the  layer  by  the  MUF  factor  M(3000).  The  term  in  parentheses  refers 
to  the  standard  ground  distance  of  3000  km  between  the  (hypothetical)  transmitter 
and  the  receiver.  The  experimental  data  for  the  numerical  factor  M(3000)  (in  terms 
of  coefficients)  come  from  13  ionosonde  stations  covering  the  geomagnetic  latitude 
range  from  71°Sto88°N. 

The  stability  and  predictability  of  the  F  layer  results  in  a  99  percent  prob¬ 
ability  (highest)  of  supporting  radio  propagation  and  communication  via  the  F  layer. 
The  next  highest  probability  is  via  the  regular  F  layer.  When  none  of  the  above 
modes  are  possible,  the  Fg  mode  is  considered  for  communication. 

For  computing  the  system  performance  the  program  allows  a  selection  from 
10  antenna  patterns.  The  program  takes  into  account  the  ground  losses,  ionos¬ 
pheric  losses,  free  space  losses,  and  the  excess  losses.  The  program  computes 
the  radio  communication  circuit  reliability,  service  reliability,  and  the  multipath 
evaluations. 

The  ITS-78  model  has  several  limitations.  The  results  from  the  model  are 
useful  only  when  the  operating  frequency  is  below  the  maximum  usable  frequency. 
The  model  assumes  that  transmission  will  be  by  reflection  from  the  ionosphere.  For 
this  the  transmitter  and  the  receiver  must  be  on  the  same  side  of  the  ionospheric 
layer  (for  example,  ground-to-ground  communication).  The  results  cannot  be 
applied  to  satellite-to-ground  communications.  The  model  does  not  take  into  ac¬ 
count  the  daytime  Fj  layer  which  usually  develops  between  the  F  and  l*^  layers. 

The  model  does  not  adequately  account  for  the  elect ron  density  above  the  altitude  of 
h  Fg-  Finally,  the  model  does  not  take  into  account  the  dependence  of  absorption 
on  the  operating  frequency  in  considering  the  D-layer  absorption. 

2.2  The  Hen!  Model 

2  14 

The  Bent  Model  #  is  basically  for  g round-to-satellite  communication  but  can 
be  adapted  for  ground-to-g round  or  satellite-to-scitellite  communication.  The  main 
purpose  of  the  Bent  model  is  to  determine  the  total  electron  content  (TFC)  of  the 
ionosphere  as  accurately  as  possible  in  order  to  obtain  high  precision  values  of  the 
delay  and  directional  changes  of  a  wave  due  to  refraction.  G round-to-satellite 
communication  demands  operating  frequencies  which  are  higher  than  the  MUF.  Thus 
the  mode  involves  the  transmission  refraction  characteristics  of  the  F^  layer  and 
the  electron  density  distribution  above  the  height  of  the  V  (}  peak  must  be  known. 

14.  Mewellvn,  S.  K.  ,  and  Bent,  R.  B.  (1973)  Documentation  and  Description  of  the 
Bent  Ionospheric  Model,  A FC RL  -  T  R- 73-0657,  AD  772733. 
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The  model  provides  (output)  the  vertical  total  electron  content  above  the  trans¬ 
mitter,  the  profile  of  vertical  electron  density  with  altitude,  and  the  total  electron 
content  along  the  path  between  the  satellite  and  the  ground.  It  also  provides  the 
refraction  corrections  to  the  elevation  angle,  the  range,  and  the  range  rate. 

The  input  parameters  to  the  model  are  the  date.  Universal  Time,  locations  of 
the  transmitter  and  receiver  (ground  and  satellite),  rate  of  change  in  elevation  and 
altitude  of  the  satellite,  operating  frequency,  the  solar  flux  (10.7  cm  flux),  and 
sunspot  number. 

The  data  base  of  the  Bent  model  consists  of  50,000  topside  ionospheric  sound¬ 
ings,  6000  satellite  measurements  of  electron  density  and  400,000  bottomside 
soundings  of  the  ionosphere.  The  data  extend  from  196  2  to  1969  to  cover  the  maxi¬ 
mum  and  the  minimum  of  the  solar  cycle. 

The  bottomside  data  are  foh^  hourly  values  from  14  stations  in  the  American 
sector  covering  geographic  latitudes  76°  N  to  12°  S  or  geomagnetic  latitudes  from 
85°  to  0°.  The  topside  soundings  cover  the  period  1962  to  1966,  with  geomagnetic 
latitude  range  85°  to  -75°,  and  the  electron  density  profiles  are  from  1000  km  down 

to  the  altitudes  of  the  foF0  peak  (h  Fn).  The  satellite  data  are  from  the  Ariel  3 

2  m2 

satellite  covering  the  period  May  1967  to  April  1968  and  are  linked  with  real  time 
foF9  from  13  ground  stations.  Thus  the  data  base  of  the  Bent  model  refers  to  solar 
cycle  20,  while  the  data  base  of  the  ITS- 7 8  model  is  from  solar  cycle  19. 

The  Bent  model  uses  foT  ^  from  the  ITS- 7 8  model.  Instead  of  the  monthly 
median  values,  the  Bent  model  computes  average  values  for  every  10-day  interval 
of  the  month  from  the  10.  7  cm  solar  flux  input.  For  the  height  h  Fg»  it  uses 
IU(3000)  factors  of  NOAA  (ITS- 78)  in  terms  of  the  sunspot  number.  It  uses  an 
empirical  polvnominal  for’  M  in  computing  h  F ' in  place  of  the  Shimazaki  equa¬ 
tion  ^  used  by  ITS- 7 8. 

The  distribution  of  electron  density  with  altitude,  assumed  bv  the  Bent  model 
for  the  computation  of  the  total  electron  density  is  shown  schematically  in  Figure  1. 
Starting  from  the  bottom,  it  divides  the  profile  into  five  sections;  a  bottom  bi¬ 
parabolic  F2  layer;  a  parabolic  F2  layer  above  the  peak;  and  three  exponential  sec¬ 
tions  to  cover  the  altitude  above  h  (h  F„  to  1000  km).  The  construction  of  the  pro- 

o  m  2  ^ 

file  needs  the  parameters  k,t  kof  k..,  Y.,  Y  ,  foFn,  and  h  Fn.  The  last  two  have 
r  1  2  3  t  m  2  m2 

already  been  explained.  The  dependence  of  the  other  parameters  on  geomagnetic 

14 

latitude,  solar  flux,  fol^*  and  the  season  is  from  work  of  Bent  .  The  topside  and 
the  first  adjoining  exponential  section  arc  matched  at  a  height  d  (above  h  )  by  the 
equation 


15.  Shimazaki,  T.  (1955)  World  wide  variations  in  the  height  of  the  maximum 

electron  density  of  the  ionospheric  Fn  layer,  J.  Radio  Res.  Labs.,  Japan 

2:86.  J 


d  ‘  h0  -  hm  '  E7  [V^V-.] 

where  Y’t  is  the  half  thickness  of  the  layer  and  k ^  is  the  exponential  constant. 


HEIGHT  h 


The  remaining  profile  above  the  F ^  peak  [of  altitude  range  (h^  +  d  to  1012)  km] 
is  divided  in  three  equal  intervals  of  altitude. 

The  model  can  predict  with  an  accuracy  of  75  -  80  percent.  If  the  model  is 
updated  with  observed  recent  data  within  a  range  of  2000  km  radius  (from  the  trans¬ 
mitters),  the  predictability  is  improved  to  90  percent. 

Though  the  model  predicts  total  electron  content  (TEC)  with  good  accuracy,  the 
model  does  not  have  separate  E  and  F^  layers.  As  the  model  was  constructed  for 
the  TEC,  the  E  and  F ^  layers  are  included  as  the  bi-parabolic  bottomside  of  the 
F 2  layer.  Also,  the  program  does  not  take  into  account  the  non-deviative  absorp¬ 
tion  in  the  underlying  D  layer. 
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2.;j  I  I i e  Ionospheric’  (loimmin  nations  Analysis  and  IVediction  Program 

(ION  CAP) 

The  lONCAP^  is  essentially  the  latest,  improved,  and  more  versatile  and 
flexible  version  of  the  ITS- 78  model. 

It  provides  80  output  options  which  can  be  divided  into  four  categories,  (1) 
ionospheric  description,  (2)  antenna  patterns,  (2)  JMLT*’  predictions,  and  (4)  system 
pe rforma nee  predictions. 

l  or  ionospheric  predictions  it  provides  monthly  median  values  for  the  param¬ 
eters  foi:,  fol  ,,  fol •  ,  h  11,  li  l‘  h  I**.,,  \  K,  \  1  ,,  and  \  l  ,,.  It  dso  pro- 

!  2  m  m  1  m2  m  ni  1  m  2 

a lesdie  er,  median,  and  upper  decile  values  of  the  minimum  of  ful  l  or  fol!  . 

It  can  also  provide  .«  prediction  in  the  form  a  plot  of  operating  frequency  with  virtual 
height  and  also  with  true  height. 

Tlie  All  1  option  of  the  output  provides  the  minimum  radiation  angle  a:ni  tile 
A!  factors  for  dl  four  modes,  K,  I  1.,,  and  K  .  i  lie  plots  for +he  diurnal  ’varia¬ 
tion  of ‘he  Ml  I  s  art  also  available.  The  Ml  Ks  provide  the  description  of  tlu 
rd.de  of  tin*  ionosphere  .»nd  do  not  include  tny  s>  stem  parameters,  i  he  operating 
i  requeue  -.  fur  a  given  radio  communication  circuit  is  the  critical  l  requeue  of  the 
1  .  rr  multiplied  b\  the  All  I  fa  do'. 

In  the  s \  '■'Ciii  pe rf< >ri na nee  options,  22  performance  parameters  are  ,c.  aiiubh. 
la-  prog."  on  :'*»*■  *.  lie  tntt  an.-  oi/put  option  computes  tin*  eh  v.dion  .agio  nd  tin 
.>:<  rditig  ncies  fo-.-  <ep*r.r.; .  dam  the  antenna  gcomet  r\  and  lh  gam. 

1  li«  •  input  :■<!•  ‘in-  program  is  4he  date,  «  im  ers.d  lime,  gei «g  r.tphimd  locations 
of  tin  *  i44:  ■  r  md  hie  r«*eei.*-r,  m<!  sunspot  iiumbe  r.  lhe  pr«*gr-.m  can  ♦vaunt 

•  \tern.l  r  pia  re.  n  i  n  <  te  m  ,s  input  to  4  he  p  rog  r  . n ; .  for  aaVnn..  pdternune 
a  M*b  >•*  4  'ii  : : * ■  nil..  *  ; *■  * • ; .ntennas  in  t:;«  program  i.  ntenn  >s  from  1  i 

bo  a  ::  ("iiiirii.  I  <  ■  r  t:  a  sVn  i  pe  rf« » r  { ;  •  nec  >ddp  a  »n  d  i  ap'd  s  S’a  r:  .  s  r..  dm  - 

t  H  -a  j.H  !  4  'a  *  . '  . !  I S  ; ! :  I  *  *  •  •  r  ,  ■  ud  *  ne  s  \  *  :t  jo  <  »f  *iie  !*e,  ei  \  e  »  re  need.  d. 

Ill-  * !  ■.  lor  M4,  e  1  e  cf  m  dc  a  -  it  \  d  i  st  ri  hut  ion  a  it  h  «H  i*  i. ue  lor  tii*  U'\c  \1’ 

ni'ogri::  i  -  >u-  r:  :u  1  ignre  ....  I  hi  model  :!  o  I  p.mabolw  !«\ers,  1!,  1.,  md  1^. 

1  la  dtd  a‘.-s  :m  *la  :>«■  .k  e  lee*  n  -a  i  ien>  jt  ies  «reii  C,  h  1, ,  md  h  t  ..  idle  half 

:u  m  1  n. 


I  .  1 

-i  1 .  11  n 


re 


•-li,  d  !  , ,  *!v!  V  J  ^  I’esneet i  .  e.l\  ,  l  or  Mu- 

;?:ad  to  he  1.  Idle  \]  J  ,\er  f  IX*  c!  a  It  it !  u  it  u 


20  Kir:.  |l  >M  A!4  i  mp  ro\  es  on  ♦  In-  MS-,  ;i  model  b\  in¬ 
i’^.  l  he  I>  i..  er  i  ont  nbut ion  m  considered  indirect^  bv 
ling  iii  o\poue:Pi  d  4  u!  Vr  the  !!  1  r.  or  down  to  the  altitude  of  70  kiiu  In  the  ‘ran 


!»* .  i  -t  i  nr  I  »  md  I  ^  \  .v. 


>ihan  region  l»  ♦  ■ .  *  *  n  the  i.  .nd  1 


i  if  I  J  IS  absent)  1  a  \  e  r  s ,  the  electron 


l‘«.  I.h, Ml,  I..,  H  colon,  I.ucas,  I).!.,  and  leters,  1  .  K.  (h'Vld 

l.stm,  Mag  die  Ib-rt'nnii.mcr  of  1  eh-comninmration  vsb-n.s  l  sing  t  he 
it im ispiie ric  Tr  ms i n i s s7<  >u  CTi  tiineV,  ST ati, mal  Tdu't > n i m u n i c ; items  and 
Triform  at  j,  m  Vdn  dnisfrd  ion~  Ri  »uh!e  r,  Colorado. 
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density  is  assumed  to  be  linear  for  the  frequency  range  f  to  f  where  f  =x  •  foE 

v  u  v  v 

and  fu  =  xu  *  foE.  Typical  values  for  and  are  0.  85  and  0.  98  respectively. 
Thus  the  electron  density  decreases  above  the  parabolic  nose  0.  85  foE  (<  foE)  and 
continues  upwards  up  to  0.  98  foE  (<foE)  producing  a  valley  in  the  transition  region. 
When  the  =  1  the  valley  is  absent  in  the  transition  region  and  the  curve  is  a 

vertical  line  starting  at  the  tip  of  the  parabolic  nose  of  the  E  layer.  The  Fj  layer 
forms  a  linear  or  parabolic  ledge  depending  on  the  magnitudes  of  h^Fj, 
foFj,  and  foF^.  In  the  ITS -7 8  model,  the  Fj  layer  is  assumed  to  be  absent.  In  the 
IONCAP  model  the  numerical  coefficients  for  foE  are  functions  of  geographic  lati¬ 
tude  for  both  solar  maximum  and  minimum  from  the  work  of  Leftin.  1  7  The  model 

1 8 

uses  foFj  maps  of  Kosich  and  Jones.  It  also  takes  into  account  the  retardation 
below  the  F^  layer. 


Figure  2.  Schematic  for 
the  Electron  Density  Profile 
and  Virtual  Height  for  the 
IONCAP  Model 


For  the  MUF  computations  the  model  uses  the  corrected  form  of  Martyn's 
theorem.  As  the  absorption  equations  using  the  secant  law  do  not  work  for  lower 


17.  Leftin,  M.  (  1976)  Numerical  Representation  of  Monthly  Median  Critical  Fre- 

quencies  of  the  Regular  E.  Region  (FnE),  OT  Report  76-88.  Moulder.  Colorado. 

18.  Rosich,  R.  K.  ,  and  Jones  W.  B.  (1973)  The  Numerical  Representation  of  the 

Critical  F requency  of  the  F  i  Region  oi  the  Ionosphere,  OT  Report  73*22. 
Boulder,  Colorado.  — 
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frequencies  at  altitudes  below  90  km,  these  equations  have  been  modified  in  the 
IONCAP  program.  The  IONCAP  provides  two  programs  (1)  the  ITS-78  short  path 
geometry  and  (2)  the  long  path  (>  10,  000  km)  geometry.  In  addition  to  the  ITS-78 
model,  the  path  computations  now  include  the  F^  mode,  the  over-the-MUF  mode, 

D  and  E  region  absorption  losses,  and  sporadic  E  losses.  A  correction  to  fre¬ 
quency  dependence  is  added  for  low  frequencies  reflected  from  altitudes  below 
<  90  km. 

The  improvements  over  ITS-78  can  be  summarized  (see  IONCAP)  as  follows: 

(1)  The  description  of  the  ionosphere  is  now  more  complete. 

(2)  The  loss  equations  have  been  supplemented.  This  includes 
E  mode  adjustments,  sporadic  E  effects,  over-the-MUF 
losses,  and  losses  for  low  reflection  heights. 

(3)  Revision  is  made  to  the  ray  path  geometry  calculations.  This 
was  an  empirical  adjustment  to  Martyn's  theorem. 

(4)  Revision  of  the  loss  statistics  was  made  to  include  the  effects 
of  the  sporadic  E  layer  and  of  over-the-MUF  modes. 

(5)  A  separate  long  path  model  is  developed. 

(6)  Revisions  were  made  to  the  antenna  gain  package. 

All  the  models  predict  only  the  quiet  ionosphere,  which  shows  a  large  system¬ 
atic  dependence  on  Latitude,  longitude,  season,  time,  and  sunspot  activity. 

2.4  The  Bradley  Model 

The  Bradley  model  contains  two  modifications  to  the  existing  models:  (1)  the 

filling  of  the  valley  between  the  E  and  F  layers  (F  ,  or  F  ifF,  is  absent),  by 
19  1  4  1 

parabolic  layers,  and  (2)  a  simple  formulation  of  the  prediction  of  the  probability 

20 

of  the  high-frequency  propagation. 

The  assumption  that  the  electron  density  distribution  in  the  E  and  F  layers  is 

parabolic  in  shape,  results  in  a  valley— a  reduction  in  electron  density  between  the 

two  peak  values,  N_E  and  N  F  — at  altitudes  h  E  and  h  F,  respectively.  The 
mm  m  m  r  J 

rocket  observations  have  shown  that  in  reality  the  electron  density  at  any  altitude 
between  the  altitudes  hmE  and  hmF  is  rarely  smaller  than  the  peak  electron  density 
N^E.  Thus  the  assumption  of  the  parabolic  shapes  for  E  and  F  layers  underesti¬ 
mates  the  electron  density  in  the  altitude  region  between  h  E  and  h  F.  To  correct 

19  m  m 

such  an  underestimation  Bradley  and  Dudeney  suggested  a  linear  distribution  of 


Bradley,  P.  A. ,  and  Dudeney,  J.  R.  (1973)  A  simple  model  of  the  vertical  dis¬ 
tribution  of  electron  concentration  in  the  ionosphere,  J.  Atmos.  Terr.  Phys. 
35:2131.  - 


20.  Bradley,  P.A.,  and  Bedford,  C.  (1976)  Prediction  of  HF  circuit  availability. 
Electronics  Letters  12:32. 
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electron  density  from  foE  to  1.  7  foE.  The  lower  end  is  at  hmE.  At  the  upper  end, 
the  F  layer  is  parabolic  in  shape  down  the  altitude  where  the  plasma  frequency  is 
1*  7  foE-  Th*s  linear  interpolation  has  not  yet  been  incorporated  in  the  IONCAP 
model  (see  Figure  2)  of  Lloyd  et  al.  16 

In  high  frequency  prediction  it  is  essential  to  know  the  probability  of  commu¬ 
nication  at  any  particular  operating  frequency  f.  For  convenience  the  observed 
data  are  expressed  as  F^— the  upper  decile,  50  percent,  and  F^— the  lower 
decile  of  the  ratios  of  f/f^.  The  distribution  functions  of  F^  and  F^  are  not  simple 
Gaussian  distributions  (Fu  -  MUF=Mui  -  F£  for  Gaussian),  The  distributions  are 
X  -distributions  of  F^  and  F^.  For  probability  determinations  these  two  x^-dis- 
tributions  (of  F^  and  F^)  have  to  be  used.  Two  variables  F^  and  F^  with  their 
associated  degrees  of  freedom,  and  the  need  to  integrate  the  x^~distribution  curve 

makes  the  process  of  determining  the  probability  distribution  very  cumbersome. 

20  u 
Bradley  and  Bedford  derived  simple  empirical  equations  for  this  probability 

distribution.  The  equations  are: 


Q 


130 


80 

T-  (f/f  ) 


l  - 


or  100,  whichever  is  smaller  for  f  ^  f 


and 


Qov  n 

=  JJJJ  p— p  ~  or  0,  whichever  is  larger  for  f  >  , 


where 

Q  -  is  the  cumulative  probability, 

f  -  is  the  operating  frequency, 

*  is  the  predicted  median  frequency, 

F c  -  is  the  lower  decile  (of  f/f  ), 

*  m 

F  -  is  the  upper  decile  (of  f/f  ). 
u  m 

They  note  that  the  probability  distribution  from  their  simple  empirical  equation 
is  as  good  as,  though  not  always  better  than,  that  from  the  \2  distribution.  There¬ 
fore  they  highly  recommend  a  replacement  of  \2  distribution  procedure  by  these 
equations  for  a  determination  of  the  present  probability  that  signals  will  propagate 
at  a  given  hour  over  a  given  sky-wave  path. 


The  latest  computer  model  like  the  TONCAP'  has  not  incorporated  'Bradley 
1 6 

Features'  in  its  program. 

2.5  The  Air  Force  Global  Weather  Central  4-D  Model 

21 

The  input  data  to  the  4-D  numerical  model  are  the  critical  frequencies  for 
the  layers  and  M(3000),  real  time  or  near  real  time  observations  from  40  ground 
stations  around  the  world,  and  total  electron  content  (TEC)  from  eleven  stations. 

The  frequency  of  observations  varies  from  hourly  (best)  at  one  end  to  weekly  (worst) 
at  the  other  end.  The  desired  purpose  of  the  4-D  model  is  to  produce  a  consistent 
ionospheric  specification  anywhere  in  the  northern  hemisphere  for  a  24-hr  period. 

In  that  sense  it  is  not  a  forecasting  model  like  the  other  models  mentioned  above. 

This  model  has  three  ionospheric  layers,  E,  F^,  and  Fg.  Each  layer  is  repre¬ 
sented  by  a  Chapman  distribution  function 


Ne<h>  s  Nemax  exp  j  at 1  -  z  -  exp  <  -  f  >1  J 


where  a  refers  to  the  loss  mechanism  and  z  is  given  by  z  =  • 


and  h  is  the 
s 


scale  height  for  the  layer.  At  a  given  altitude  the  total  contribution  to  electron 
density  is  the  sum  of  contributions  by  all  the  three  layers. 

For  any  height  the  electron  density  is  approximated  by 
N 

N  <f)  =  E  a.W.  (t) 
k=  1 

where  a  is  the  weighting  factor  and  W(f )  is  an  empirically  derived  set  of  discrete 
orthogonal  functions  for  the  altitude  interval  ’ f  ' .  The  95  to  2000  km  range  is 
divided  into  127  intervals.  The  widths  of  the  intervals  range  from  5  km  at  the  lowest 
altitudes  to  50  km  at  highest  altitudes.  The  empirically  derived  function  W(£)  is  in 
two  parts,  spherical  harmonic  functions  for  spatial  dependence  and  trigonometric 
functions  (sine,  cosine)  for  temporal  dependence. 

With  the  help  of  these  variables  and  W^,  the  entire  data  base  for  the  ionos¬ 
phere  is  reduced  to  a  limited  number  of  coefficients.  These  can  be  used  to  construct 
the  electron  density  profile,  for  any  location  in  the  Northern  Hemisphere,  valid  for 
a  24-hr  specification  period.  The  model  is  still  being  developed.  The  specification 
accuracy  of  the  model  will  strongly  depend  on  the  frequency  and  reliability  of  the 
input  data —real -time  experimental  observations  from  the  40  ground  stations  of  the 

21.  Flattery,  T.W.,  Tascione,  T.F.,  Secan,  J.A.  ,  and  Taylor,  Jr.,  J.  W. 

A  Four-Dimensional  Ionospheric  Model  (private  communication). 
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northern  hemisphere.  Also  the  quality  of  specifications  interpolated  for  locations 
inside  the  network  will  be  better  than  those  extrapolated  outside  the  observational 
network. 

3.  THE  THEORETICAL  MODELS 

The  theoretical  models  for  the  ionosphere  are  based  on  the  physical  processes 
responsible  for  the  observed  ionospheric  phenomena.  The  processes  responsible 
for  the  ionosphere  are  production,  maintenance,  and  decay  of  the  ionosphere.  As 
the  approach  deals  directly  with  the  physical  processes,  and  not  with  the  observed 
phenomena,  the  emerging  model  is  called  a  physical  model. 

Four  models  are  summarized  in  Table  1  to  show  several  variations  in  the  same 
processes  considered  by  different  workers.  Strobel  and  IVlcElroy9  considered  only 
the  F2  region  (200  to  7  00  km),  whereas  others  took  into  account  the  altitude  range 
from  120  to  1200  km.  Nisbet10  constructed  the  first  computer-based  simple  physical 

model  JV1K-1  for  the  ionosphere.  He  considered  only  three  neutral  constituents  Nn, 

8  12  2 
C>2  and  O,  whereas  Stubbe  and  Oran  and  Young  also  considered  the  minor  con¬ 
stituents  He  and  H.  For  the  dissociation  and  ionization  of  the  neutral  species,  the 
incident  solar  EUV  radiation  in  the  range  30  to  1912  A  is  used,  along  with  the 
wavelength  dependent  absorption  and  ionization  cross  sections  for  the  neutral 
species.  Nisbet10  considered  three  basic  predominant  ionic  species:  G+,  NO+,  and 
O/.  Oran  and  Young12  took  into  account  the  additional  ionic  species  N9+,  He+,  N+, 
and  H  .  One  has  to  consider  the  chemical  reactions  which  produce  ions  by  charge 
exchange  processes.  Nisbet19  used  5  reactions  whereas  Oran  and  Young12  used 
24  chemical  reactions  (see  Ref.  9).  For  maintenance  of  the  ionosphere  the  pro¬ 
cesses  of  diffusion  and  photo-ionization  are  assumed.  The  processes  of  dissociative 
recombination  and  radiative  recombination  are  responsible  for  the  decay  of  the 
ionosphere.  For  his  simple  model,  Nisbet  neglected  the  transport  processes  such 
as  neutral  winds,  electric  fields,  and  magnetic  fields.  The  procedure  is  further 
c  omplicated  because  coupled  simultaneous  equations  must  be  solved  for  neutral 
winds,  mass  transport,  and  energy  transport.  For  determining  electron  density  in 
the  ionosphere,  the  gas,  consisting  of  both  ions  and  electrons,  is  considered  elec¬ 
trically  neutral.  Thus,  in  every  elementary  volume,  the  number  of  electrons  is 
equal  to  the  number  of  ions.  All  the  models  reproduce  many  of  the  observed 
features  such  as  the  diurnal  variation,  seasonal  variation,  and  solar  cycle  depend¬ 
ence  of  the  mid-latitude  ionosphere  under  quiet  conditions.  The  accuracy  of  the 
theoretical  models  depends  upon  the  understanding  of  the  physical  processes  con¬ 
sidered  in  the  models.  For  accurate  predictions  from  the  theoretical  models,  pre¬ 
cise  information  on  the  large  number  of  variables  used  in  the  models  is  necessary. 


Also,  the  models  use  several  observed  average  boundary  conditions  which  could 
have  a  large  variability  dependent  on  other  geophysical  parameters  such  as  solar 
activity  and  magnetic  activity.  The  results  from  the  models  are  adequate  for  long 
term  planning  of  science  and  engineering  applications.  Though  these  models  re¬ 
produce  main  observed  average  features  of  the  ionosphere,  these  are  unable  to 
specify  the  ionosphere  within  an  accuracy  of  ±  20  percent  needed  by  the  systems  in 
operation.  At  present,  the  main  input  information  of  solar  EUV  radiation  needed 
for  the  theoretical  models  is  not  routinely  available  for  predicting  the  ionosphere. 


4.  COMPARISON  OF  THE  PHENOMENOLOGICAL  MODELS, 

THEIR  LIMITATIONS  AND  AVAILABILITY 

In  comparing  the  models  one  must  note  that  'IONCAP'  is  the  modified  and  more 
flexible  version  which  replaces  the  ITS-78  model.  As  the  ITS-78/IONCAP  and  the 
Bent  models  serve  entirely  different  purposes,  it  is  essential  to  understand  the 
difference  in  their  approaches  and  final  output  parameters  computed  by  the  models. 
These  are  summarized  in  Table  2.  The  left-hand  column  in  Table  2  lists  the 
parameter  under  consideration.  The  next  three  columns  summarize  the  features 
in  each  of  the  models,  ITS -7 8,  IONCAP,  and  the  Bent  model  respectively.  From 
the  table  it  is  seen  that  the  selection  of  a  model  will  depend  more  upon  the  informa¬ 
tion  sought  under  the  parameter  headings  than  on  accuracy.  The  IONCAP  model  is 
basically  useful  for  wave  propagation  using  operating  frequencies  which  would  be 
reflected  by  the  E,  Eg,  F^,  and  layers.  On  the  other  hand,  the  Bent  model 
relies  on  the  transmission,  refraction,  and  absorption  characteristics  of  the  ionos¬ 
phere,  with  the  operating  frequency  much  larger  than  the  foFg  frequency.  All  the 
models  predict  quiet  ionospheric  conditions  only.  The  models  do  not  hold  for  dis¬ 
turbed  ionospheric  conditions. 

The  additional  limitations  of  these  models  are: 

(1)  Ail  the  models  are  poor  in  predicting  the  high  latitude  ionosphere, 

(2)  None  of  the  models  take  into  account  the  effects  of  particle 
precipitation  in  the  auroral  region  which  enhance  the  E  (Eg) 
and  F  layers, 

(3)  The  mid-latitude  trough  which  exhibits  large  horizontal 
gradients  in  electron  density  is  not  incorporated  in  these 
models, 

(4)  These  models  are  good  for  latitudes  ±  20  to  ±  60°,  and  are 
poor  predictors  for  the  equatorial  region  and  the  high  latitude 
region. 
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Nonetheless  these  models  serve  two  useful  functions:  (1)  to  predict  ionospheric 
parameters;  and  (2)  to  determine  physical  phenomena  and/or  to  modify  existing 
coefficients  for  explaining  the  deviations  between  the  experimentally  observed  value 
and  the  predictions  from  these  models. 

1  16 

The  computer  programs  for  the  ITS-78  and  the  IONCAP  program  are  avail¬ 
able  from  the  Institute  for  Telecommunications  Sciences,  Boulder  Colorado  80303. 

14 

The  computer  programs  for  the  Bent  model  are  available  from  the  Atlantic 
Science  Corporation,  P.  O.  Box  3201,  Indialantic,  Florida  32903. 


Table  2.  Intercomparison  of  the  Empirical-Computer  Based  Ionospheric  Models 


Parameter 

Ionospheric  Models 

ITS -7  8 

IONCAP 

Bent 

D-Region 

Non -deviative  and 
deviative  absorp¬ 
tions  only 

Same  as  ITS  -7  8 
+  E  Layer 
exponential 
extension  down 
to  7  0  km 

Not  modeled 

E  -Region 

foE 

hmE 

YmE 

Modeled  by 

Leftin  et  al^ 

110  km  Fixed 

20  km 

Parabolic  Shape 

Same  as  ITS -7  8, 

+  exponential 
down  to  70  km. 
Leftin^7 
coefficients 

Not  modeled 

F ^  -  Region 
foF  j 

Not  modeled 

Rosich  &  Jones 
coefficients 

YrnF  j 

m  1 

h  F  .  /  F  ,  -  4 

m  1  ym  1 

(fixed) 

Not  modeled 

F2-Region  Bottomside 

foF2 

hmF2 

Y  F 
m  2 

Hayden  Lucas 
coefficients^ 
Shimazaki  eq, 

+  E  layer 
retardation 

Kelso^ 

Same  as  ITS-78 

Bi  -parabolic 

Bent  coefficients. 

F2 -Region  Topside 

Not  modeled 

Not  modeled 

Up  to  1000  km 

(Due  to  large  number  of  references  cited  in  Table  2,  they  will  not  be  listed  here. 
See  References,  page  25.  ) 


Table  2.  Intercomparison  of  the  Empirical -Computer  Based  Ionospheric  Models 
(Cont) 


Parameter 

Ionospheric  Models 

ITS -7  8 

IONCAP 

Bent 

E-F  Transition  Region 

Not  modeled 

Modeled 

Not  Modeled 

Electron-Density  Profile 

Not  computed 

Available  up 

to  h  FfJ 
m  2 

Available  up  to 
3500  km 

Total  Electron  Content 
(TEC) 

Not  computed 

Not  computed 

Com  puted 

MUF 

For  short  path 
only 

Also  for  long 
path 

<>  10,000  km) 

Not  modeled 

Short-Term  Prediction 
of  MUF 

Function  of  Kp 

Not  modeled 

Not  modeled 

Input  Parameters 
required 

Sunspot  number 

Sunspot  number 

Sunspot  number 
and  10.  7  cm 
solar  flux 

Noise  parameters 

Galactic 

Same 

Not  modeled 

Atmospheric 

Manmade 

Modified 

Same 

MUF  50% 

EOT  90%  Fu  * 

HPF  10%  F^^ 

Modeled 

Modeled 

Not  modeled 

System  Performance 

Modeled  for 
short  path 
(<  3000  km) 

Also  has  a  long 
path 

option  - 
10,  000  km 

Not  modeled 

Antenna  Patterns 

Uses  ITSA-1 

Package  with 

10  antenna 
options 

Modified 

ITS-78 
package  with 

17  antenna 
options 

Not  modeled 

Sporadic  E 

Modeled  in  terms 
of  occurrence 
frequency 

Same  as  ITS-78 

Not  modeled 

Circuit  Reliability 

Service  Probability 
Multipath  Evaluation 

Modeled 

Modeled  as 

ITS-78 

Not  modeled 
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L  LOSS  A  RY 


foL 
foF  j 
foF9 
foK 

s 

f 

Ff 

fm 

F 

u 

FOX 

h 

h 

s 

h 

max 

h  E 
m 

h  F 
m  1 

h  F0 
m  2 

k.,  k.,,  k.. 


LUF 


MU' 


N  (h) 
e 

N 

cmax 


N 

m 


E 


N  F 


m 

vvk(n 


Y 

Y 

Y 


m 


E 


m 

m 


h  2 1 


V 


t 


Weighting  factor  in  expression  for  electron  density  as 
a  function  of  height. 

Loss  coefficient  in  the  Chapman  electron  density  profile. 
E-layer  critic.il  frequency. 

F  ^ -  layer  critical  frequency. 

F9 -layer  critical  frequency. 

Sporadic  E-laver  critical  frequency. 

Operating  frequency  of  communications  system. 

Lower  decile  of  f/f  . 

Median  value  of  Maximum  Usable  Frequency. 

Upper  decile  of  f/f  . 

m 

Optimum  frequency  for  transmission. 

Height. 

Scale  height  of  layer. 

Height  of  maximum  electron  density. 

Height  of  maximum  electron  density  within  the  E-layer. 
Height  of  maximum  F^ -layer  electron  density. 

Height  of  maximum  F^-layer  electron  density. 
Exponential  constants  in  the  Bent  model  for  the  three 
portions  of  the  electron  density  profile  above  fol'^. 
Lowest  usable  frequency. 

Maximum  usable  frequency. 

Electron  density  as  a  function  of  height. 

Maximum  electron  density. 

Peak  E-layer  electron  density. 

Peak  F- layer  electron  density. 

Set  of  empirically  derived  discrete  orthogonal  functions 
for  the  altitude  interval  C  in  the  AFGWC  4-D  model. 

Half  thickness  of  the  E  layer. 

Half  thickness  of  the  Fj  layer. 

Half  thickness  of  the  F ^  layer. 
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